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Abstract 

The prediction of thermo-mechanical behaviour of heterogeneous materials 
such as heat and moisture transport is strongly influenced by the uncer- 
tainty in parameters. Such materials occur e.g. in historic buildings, and the 
durability assessment of these therefore needs a reliable and probabilistic sim- 
ulation of transport processes, which is related to the suitable identification 
^ ! of material parameters. In order to include expert knowledge as well as ex- 

Ti^ij- \ perimental results, one can employ an updating procedure such as Bayesian 

O ■ inference. The classical probabilistic setting of the identification process in 

Bayes's form requires the solution of a stochastic forward problem via com- 
putationally expensive sampling techniques, which makes the method almost 
impractical. 

In this paper novel stochastic computational techniques such as the stochas- 
tic Galerkin method are applied in order to accelerate the updating proce- 
^ ■ dure. The idea is to replace the computationally expensive forward simula- 

. tion via the conventional finite element (FE) method by the evaluation of a 

polynomial chaos expansion (PCE). Such an approximation of the FE model 
for the forward simulation perfectly suits the Bayesian updating. 

The presented uncertainty updating techniques are applied to the numer- 
ical model of coupled heat and moisture transport in heterogeneous materials 
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with spatially varying coefficients defined by random fields. 
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1. Introduction 

Durability of structures is infiuenced by moisture damage processes. High 
moisture levels cause metal corrosion, wood decay and other structural degra- 
dation. Thermal expansion and contraction, on the other hand, can induce 
large displacements and extensive damage to structural materials with dif- 
fering coefficients ,e.g. masonry. The Charles Bridge in Prague, currently 
the subject of rehabilitation works, is a typical example, see jH. A study 
of the coupled heat and moisture transport behaviour is thus essential in 
order to improve the building materials' performance. So far, a vast number 
of models have been introduced for the description of transport phenomena 
in porous media. An extensive overview of transport models can be found 
in Q . In this work we focus on the model by Kiinzel jsj , since the predicted 
results comply well with the results of experimental measurements j^, once 
the relevant material parameters are well estimated. 

Material properties are usually determined from experimental measure- 
ments via an identification procedure, see e.g. jsj. However, the experimental 
measurements as well as the identification methods involve some inevitable 
errors. Bayesian updating, employed within this study, provides a general 
framework for inference from noisy and limited data. It enables mutually in- 
volving both expert knowledge of the material, such as limit values of physical 
parameters, and information from experimental observations and measure- 
ments. In other words, it uses experimental data to update the so-called 
a priori uncertainty in the material description and results in a posterior 
probabilistic description of material performance 0. In addition, unlike tra- 
ditional identification techniques that aim to regularise the ill-posed inverse 
problem to achieve a point estimate, the Bayesian identification process leads 
to a well-posed problem in an expanded stochastic space. 

The main disadvantage of Bayesian updating lies in the significant com- 
putational effort that results from the sampling-based estimation of posterior 
densities |7|]. While deterministic quadrature or cubature may be attractive 
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alternatives to Monte Carlo at low to moderate dimensions [8|, computa- 
tionally exhaustive Markov chain Monte Carlo (MCMC) remains the most 
general and flexible method for complex and high-dimensional distributions 



[9|, 110[. In a sampling-based procedure, the posterior distribution must be 
evaluated for any sample generated from the prior one in order to decide, 
whether the sample is admissible or not. The computation of the posterior in- 
volves the evaluation of the computational model — the FE discretisation of a 
non-linear partial differential equation (PDE) — relating model (i.e. material) 
parameters and observable quantities (i.e. model outputs). Hence, complex 
and time-consuming models can make the sampling procedure practically 
unfeasible. 

Bayesian updating of uncertainty in the description of the parameters of 
Kiinzel's model is thoroughly described in ll| for the case of heterogeneous 
material, where material parameters are described by random flelds (RFs). 
It was shown that Bayesian updating is applicable even for such a complex 
and nonlinear model as Kiinzel's model. However, the demonstrated results 
were performed for a sample with a coarse FE, thereby rendering the eval- 
uation of the numerical model computationally relatively cheap. A higher 
complexity of modelled structure and its FE-based numerical model lead to 
time-consuming simulations and are prohibitive for the sampling procedure. 
In such a case one may construct an approximation of the model response and 
evaluate this within the sampling procedure in order to render the updating 
procedure feasible [l2|, [l3 . 

The efficient forward propagation of uncertainty, which may describe ma- 
terial properties, the geometry of the domain, external loading etc., from 
model parameters to model outputs is a main topic of stochastic mechanics. 
The recently developed polynomial chaos (PC) variant of the stochastic finite 
element method (SFEM)— the spectral SFEM (SSFEM) Hi, [lil, [l6|, [13, [l8 — 



has become one of the promising techniques in this area. Some of the uncer- 
tainties in the model are represented as random fields/processes. Here one 
often employed technique in SFEM computations is the use of a truncated 
Karhunen-Loeve expansion (KLE) to represent the RFs in a computation- 
ally efficient manner by means of a minimal set of random variables (RVs) 



19L Il6l . |20| . |21[, via an eigenvalue decomposition of the covariance. This 
approach involves the introduction of an orthogonal — hence uncorrelated — 
basis in a space of RVs. These are projections of the RF onto the orthogonal 
KL eigenfunctions, and in the case of Gaussian RFs consists of Gaussian 
RVs. In that case they are not only uncorrelated but independent — a com- 
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putationally very important property [22*]. However, the material properties 
very often cannot be modelled as Gaussian due to crucial constraints such as 
positive definiteness, boundedness in some interval, etc. In such a case, one 
has to adopt non-Gaussian models and their corresponding approximations, 
see 171, ll8|, |23|, often as a non-linear transformation of a Gaussian RF. The 
orthogonal or uncorrelated RVs alluded to above are not Gaussian in that 
case, and hence not independent. One then may adopt a pure PC repre- 
sentation of the RF in terms of polynomials of independent Gaussian RVs 



14j |. or — to take advantage of the dimension reduction inherent in the KLE 



truncation — one uses the PC representation for the orthogonal/uncorrelated 



non- Gaussian RVs from the KLE 16, 20 



In this paper, we focus on Kiinzel's model [2J, |3|, defined by uncertain 



positive-definite material parameters, modelled as log-normal RFs according 
to the maximum entropy principle. Since these RFs are non-Gaussian, their 
spectral decomposition (KLE) gives a set of uncorrelated but not necessarily 
independent RVs. To address this problem, we project the RVs onto a PC 
basis constructed from Hermite polynomials in independent Gaussian RVs as 
alluded to in the previous paragraph. Such a combined expansion (KL/PC) 
is then used to represent the RFs as inputs to the FE discretisation of the 
nonlinear Kiinzel model. The solution procedure of Galerkin type for this 
SPDE is chosen in an "intrusive" manner based on analytic computations in 



the PC/Hermite algebra [20|, |25|, |26[ . This brings huge computational savings 
in case of small and moderate problem dimensions, but it requires complete 
knowledge of the model (the FEM system can not be used in black-box 
fashion) . 

Once such a representation is propagated through the physical model, 
one obtains a description of all desired output quantities in terms of simply 
evaluable functions — in this case polynomials — of known independent Gaus- 
sian RVs. This is often called a surrogate model or a response surface. 

The paper is organised in the following way. The next Section |2] reviews 
Kiinzel's model. Section [3] is focused on the probabilistic description of het- 
erogeneous material properties where particular material parameters are not 
spatially constant. Intrusive stochastic Galerkin method for computing coef- 
ficients of the PC-based surrogate of outputs of Kiinzel's model is developed 
in Section H] and the related outcomes a presented in Section [51 Finally, Sec- 
tion |6] presents the Bayesian updating procedure on Kiinzel's model with the 
results summarized in Section [71 and Section [81 concludes. 
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2. Coupled heat and moisture transfer 



Kiinzel |3| derived balance equations describing coupled heat and 



moisture transport through porous media using the concepts of Krischer 



and Kiessl. Krischer [27[ identified two transport mechanisms for material 
moisture, one being the vapour diffusion and the other being described as 
capillary water movement. In other words, he introduced the gradient of 
partial pressure in air as a driving force for the water vapour transport and 
the gradient of liquid moisture content as the driving force for the water 



transport. This model is then extended by Kiessl [28[ who introduced the so- 
called moisture potential $ used for unification of the description of moisture 
transport in the hygroscopic < 0.9 and over-hygroscopic if > 0.9 range 
(where ip is relative humidity). The introduction of the moisture potential 
brings several advantages, especially very simple expressions for the moisture 
transport across the interface. On the other hand, the definition of the 
moisture potential in the over-hygroscopic range was too artificial, and Kiessl 
introduced it without any theoretical background, see 

For the description of simultaneous water and water vapour transport 
Kiinzel chose the relative humidity (f as the only moisture potential for both 
the hygroscopic and the over-hygroscopic range. He also divided the over- 
hygroscopic region into two sub-ranges — the capillary water region and su- 
persaturated region — where different conditions for water and water vapour 
transport are considered. In comparison with Kiessl's or Krischer's model 
Kiinzel's model brings certain simplifications. Nevertheless, the proposed 
model describes all substantial phenomena and the predicted results comply 
well with experimentally obtained data j^. Therefore, it was chosen as a 
physical basis for the formulation of the probabilistic framework. 

Kiinzel's model is described by the energy balance equation 

H M r\f) 

— - = v^[A(¥.) v^] + K{e)v^Me)v{^p..,m (i) 

and the conservation of mass equation 

= + V^[5p(0)V{mat(^^)}] , (2) 

where the transport coefficients defining the material behaviour are nonlinear 
functions of structural responses, i.e. the temperature 9[°C] and moisture 
fields. We briefiy recall their particular relations [s]: 
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Thermal conductivity [Wm K 



A = Ao|l + ''-'"'"'-'''^ ') . (3) 

Ps{b - 



Evaporation enthalpy of water [Jkg ^]: 



273 15 ^ (0.267+3.67- lo-^'e) 



• Water vapour permeability [kgm~^s^^Pa^^]: 

= (9 + 273.15)-. (6) 

• Water vapour saturation pressure [Pa]: 

• Liquid conduction coefficient [kgm^^s^^]: 

n = 3.8— . io(f-^)(-f-i) . -J^ -A . (7) 

Wf (o — v?)^ 

• Total enthalpy of building material [Jm~^] : 

H = p,c,e. (8) 

• Water content [kgm~^]: 

w = Wf—^ — . (9) 

b — (f 

A more detailed discussion on the transport coefficients can be found in 
js], lioj. Some of them defined by Eqs. ©-([H]) depend on a subset of the 
material parameters listed in Tab. [H The approximation factor b appearing 
in Eqs. ([3]) and can be determined from the relation: 

^^O^S^Wso^^ (10) 
Wso — 0.8w{ 
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where Wso is the equihbrium water content at 0.8 [— ] relative humidity. More- 
over, the free water saturation Wf must always be greater than wsq. Therefore 
we introduce the water content increment dwf > and define the free water 
saturation as 

W{ = wso + dwf . (11) 

Consequently, Wsq and dw{ substitute b and W{ as material parameters to be 
identified within the updating procedure. Tab. [T] presents the resulting list 
of W = 8 material parameters to be identified. As an outcome of such a 
substitution, all identified parameters should be positive and thus described 
by log-normal RFs (a priori information) with second order statistics (mean 
values /ig and standard deviations aq) given in Tab. [TJ Those particular 



values are chosen to correspond to materials used in masonry [30 



Parameter 






Pq 


^9 


dwf 


[kgm^ 




water content increment 


100 


20 


wso 


[kgm^ 




water content at 0.8 [— ] relative humidity 


50 


10 


Ao 


[Wm- 




thermal conductivity of dry material 


0.3 


0.1 


^tcs 


[-] 




thermal conductivity supplement 


10 


2 




[-] 




water vapour diffusion resistance factor 


12 


5 


a 


[kgm" 


-2g-0.5] 


water absorption coefficient 


0.6 


0.2 


Cs 


[Jkg- 




specific heat capacity 


900 


100 


Ps 


[kgm" 




bulk density of building material 


1650 


50 



Table 1 : Mean values and standard deviations of material parameters 

The partial differential equations and ([2]) are discretised in space by 
standard finite elements. This also goes well with the use of the stochastic 
Galerkin method for the discretisation in the stochastic space. Performing 
first only the spatial discretisation, the temperature and moisture fields are 
spatially approximated as 

N N 

9{x) = ^(f)nix)ue,n, ^{x) = ^0„(a;)M^,n (12) 

n=l n=l 

where N is the number of nodes in FE discretisation, (f)n{x) are the shape 
functions (according to the type of used elements) and Me,n and M(^,„ are the 
nodal values of temperature field 9 and moisture field ip, respectively. 
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Using the approximations Eq. ( IT^ and Eqs. ([T]), ([2]), we obtain a set of 
first order differential equations 

fill 

K{u)u + C{u)^ = F, (13) 

where K{u) is the conductivity matrix, C{u) is the capacity matrix, = 
{ug^i, . . . , ug^N, u<p,i, • • • , M(^,iv) is the vector of nodal values, and F is the vector 
of prescribed fluxes transformed into nodes. For a detailed formulation of the 
matrices K{u) and Cju) and the vector F, we refer the interested reader 
to the doctoral thesis [Si!, Chapter 3.1]. 

The numerical solution of the system Eq. (1131) is based on a simple tem- 
poral finite difference discretisation. If we use time steps Ar and denote the 
quantities at time step i with a corresponding superscript, the time-stepping 
equation is 

n^+i = n^ + Ar[(l-7)n* + 7tt^+^], (14) 

where 7 is a generalised midpoint integration rule parameter. In the results 
presented in this paper the Crank-Nicolson (trapezoidal rule) integration 
scheme with 7 = 0.5 was used. Expressing vJ'^^ from Eq. f|T^ and substitut- 
ing into the Eq. f|T3|) . one obtains a system of non-linear equations: 

(7ArK*+i + C'+^)u'+^ = 7ArF*+i + C'^^[u' + Ar(l - -f)u% (15) 

which can be solved by some iterative method such as Newton-Raphson. 
For clarification and easier reading, we rewrite Eq. f lTSj) using the symbols 
^AtK'+\u'+^) + C'+\u'+^) and f+\u'+^) := -fArF'^^ + 
C*+^(w*^^)[w' + Ar(l — 7)1*'] in the following form 

A'+\u'+')u'+' = f+\u'+^). (16) 

3. Uncertain properties of heterogeneous materials 

When dealing with heterogeneous material, some material parameters 
can vary spatially in an uncertain fashion and therefore RFs are suitable for 
their description. This means that the uncertainty in a particular material 
parameter q is modelled by defining q{x) for each a; G ^ as a RV q{x) : i? — )■ 
M on a suitable probability space (J7,^,P) in some bounded admissible 
region ^ C M'^. As a consequence, g:^xi7— >]Risa RF and one may 
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identify i? with the set of all possible realisations of q. Alternatively, q{x, uj) 
can be seen as a collection of real- valued RVs indexed hy x E Q. 

The description of log-normal RFs given in Tab. [T]can be derived from a 
Gaussian RF g{x,u), which is defined by its mean 

figix) = E[g{x, u)] = [ gix, u) ¥{du) , (17) 

and its covariance 

Cg{x,x') = E[{g{x,u) - fig{x)){g{x',u) - fig{x'))] 

{gix, cu) - fig{x))ig{x', u) - fig{x')) P(da;) . (18) 



The log-normal RF q{x, u) can be then obtained by a nonlinear transforma- 



tion of a zero-mean unit-variance Gaussian RF g{x, u) [20|, |26 



as 



q{x,u) = exp{ij,g{x) + agg{x,uj)) . (19) 

The statistical moments Hg and cXg of the Gaussian field can be obtained 
from the statistical moments fig and ag given for the log-normally distributed 



material property according to the following relations [26 




1 



In numerical computation random fields are first spatially discretised by 
finite element method (see Eq. (fT2l) ) into a finite collection of points {a:;"^^} e 
Q. Further, the semi-discretised RF are described by a finite — but probably 
very large — number of RVs q{uj) = {q{xi, u), . . . , q{Xn, u)), which are usually 
highly correlated. Large number of RVs is, however, very challenging for the 
efficient numerical implementation of forward problem, as well as for MCMC 
identification. As already alluded to previously, the number of RVs can 
be reduced by the approximation g{uj) of a RF gjoj ) based on a truncated 



KLE including much smaller number of RVs [20|, ll3| . Here we use the KLE 
on the underlying Gaussian field g{uj), and hence the RVs in the KLE are 
independent Gaussian RVs, as already indicated above. 

The spatial discretisation of a given RF concerns also the discretisation 
of corresponding covariance function Cg{x, x') into the covariance matrix Cg 
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which is symmetric and positive definite |16l . |20| . The KLE is based on the 
spectral decomposition of the covariance matrix Cg leading to the solution 
of a symmetric matrix eigenvalue problem 



C,V. = , (21) 

where il^i are orthogonal eigenvectors and are positive eigenvalues ordered 
in a descending order. The KLE approximation g{uj) of a RF g{u:) can then 
be written as 

M 

g{uj) = /Xg + ^ Q '^i(w)V'i, (22) 

i=0 

where = tl^i'^^giu) — ^ig)/<,i are uncorrelated RVs of zero mean and unit 
variance, and in case that g{x, uj) and hence g{uj) are Gaussian, then Xiiu) are 
Gaussian and independent. The number M < n — the number of points used 
for the discretisation of the spatial domain — is chosen such that Eq. (1221) gives 
a good approximation, i.e. captures a high proportion of the total variance. 
Higher values of M lead to better description of a RF, smaller values imply 
faster exploration by MCMC. The eigenvalue problem Eq. fl2T]) is usually 
solved by a Krylov subspace method with matrix approximation. 

For large eigenvalue problems, the authors in |33] propose efficient low-rank 
and data sparse hierarchical matrix techniques. The approximation of a 
non-Gaussian RF can be then obtained by a nonlinear transformation of 
the KLE obtained for a Gaussian RF such as in our particular case, where 
the approximation of a given RF q{oj) is obtained from the Eq. (fT9|) by the 
substitution of the Gaussian RF g{uj) by its KLE giu). 

We assume full spatial correlation among material properties, i.e. spatial 
fiuctuations for all parameters differ only in magnitude. Taking into account 
a log-normal distribution of the parameters, the final formulation of the RF 
describing the parameter q then becomes 

q{u) = exp ^fig + ^3 5^ yQ6(^)V'i^ , (23) 

where the exponential is to be used at each spatial point, i.e. for each com- 
ponent of the vector inside the parentheses. The statistical moments fig and 
ag are derived from the prior mean fiq and standard deviation aq for each 
material parameter according to Eq. (1201) . The eigenvectors ipi are obtained 
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for the a priori exponential covariance function 



C{x, x') = exp 



r2| 



''X2 



(24) 



where r = (ri,r2) = x — x', and Ix^ = 0.1 [m] and — 0-04 [m] are a 
priori covariance lengths. Determination of correlation lengths is generally 
not obvious. In material modelling, one possible way is based on image 
analysis as described in 33| . A numerical study for a differing number of 



modes M included in the KLE is presented in [11 



4. Surrogate of Kiinzel's model 

While the KLE can be efficiently applied to reduce the number of RVs 
and thus to accelerate the exploration of the MCMC method in terms of 
the number of samples, construction of a surrogate of the computational 
model can be used for a significant acceleration of each sample evaluation. 



In |12l . |13[ methods were introduced for accelerating Bayesian inference in 
this context through the use of stochastic spectral methods to propagate 
the prior uncertainty through the forward problem. Here we employ the 



stochastic Galerkin method |15l . |16| to construct the surrogate of Kiinzel's 
model based on polynomial chaos expansion (PCE). 

According to Eq. ( 123|) . all model parameters are characterised by M in- 
dependent standard Gaussian RVs ^{u) = [C^i^u), . . . ,^m{^)]- Hence, the 
discretised model response u{^{lo)) = (. . . ,Ui{^{u)) is a random vec- 
tor which can be expressed in terms of the same RVs ^(w). Since ${u}) are 
independent standard Gaussian RVs, Wiener's PCE based on multivariate 
Hermite polynomials — orthogonal in the Gaussian measure — {Ha{^{uj))}aej 
see I20I for the notation) is the most suitable choice for the approxima- 



tion ti(^(u;)) of the model response ^^(^(a;)) |3J], and it can be written as 

ui^iu)) = J2^M^{co)) (25) 

where Ua is a vector of PC coefficients and the index set J' C Nq^"* is a 
finite set of non-negative integer sequences with only finitely many non-zero 
terms, i.e. multi-indices, with cardinality \J'\ = R. We collect all the PC 
coefficients in u := [. . . , u^, ■ ■ ■ ]a&j- Assuming the uncertainty in all material 
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parameters listed in Tab. [T]and consequently in the model response, Eq. f ll6p 
can be rewritten as 



A^+i u^^\^)) u^^\^) = u^^\^)) . (26) 

Substituting the model response by its PC approximation ii^'^^i^) 

given in Eq. fl25|) and applying a Bubnov-Galerkin projection, one requires 
that the weighted residuals vanish: 

W^ej: E([r+i - A'^\^)u'^\m(^m = , (27) 

where A*"^\^) := A'+^ n*+^(|)). Eq. ([27]) together with Eq. leads to 
V/3 G J : (/7,(0A'^^'(^)ifJ0) <+i = Eif^'m^m , (28) 

which is a non-linear system of equations of size N x R. 

The approximation 'U*"'"^(^) can be represented through its PC coeffi- 
cients u*"*"^, and similarly for all other quantities. Denoting the block-matrix 
A^+i(u^+i) := (E(i7;3(^)A^+^(0^«(0))/3,aej, and the right hand side P+i := 
(E(/'+^(^)i7^(^)))^gj-, the system Eq. (125]) may succinctly be written as 

A*+i(u^+i)u'+i = r+\ (29) 

The matrix A*"*"^ has more structure than is displayed here, but this is outside 



the scope of this paper; see [IG], l20| for details and possible computational 
procedures. 

The evaluation of expected values in Eq. fl28]) can often be performed 
analytically in intrusive Galerkin procedures — that is their advantage — using 
the Hermite algebra In case they are to be computed numerically, they 
may be approximated by a weighted sum of samples drawn from the prior 
distributions. To that purpose, one can apply some integration technique: 
the Monte Carlo (MC) method, the quasi-Monte Carlo (QMC) method, or 



some quadrature rule, see [20| for a recent review. The latter ones allow to 
take advantage of a possibly regular behaviour in the stochastic variables and 
consequently reduce the number of samples. Since the system of equations 
Eq. fl28]) can be quite large, the evaluation of the left hand side for each sample 
of ^ becomes costly. Here we apply a sparse-grid Smolyak quadrature rule 
sometimes also named hyperbolic cross integration method, 
which is an efficient alternative for integration over Gaussian RVs. 
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After solving the system Eq. fl29p , one has via Eq. (123]) a surrogate repre- 
sentation of the model outputs. This model approximation may be evaluated 
orders of magnitude more quickly than the evaluation containing the full FE 
simulation. 



5. Numerical results for the uncertainty propagation 

For an illustration of the described method, we employ the same simple 
example as in ll| with the two-dimensional rectangular domain discretised 
by an FE mesh into N = 80 nodes and 120 triangular elements. Its geometry 
together with the specific loading conditions are shown in Fig. [H The initial 
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Figure 1: Experimental setup 



temperature is 9in = 14 [°C] and the moisture ipin = 0.5 [-] in the whole 
domain. One side of the domain is submitted to exterior loading conditions 
^cxt = 5 [°C] and ipext = 0.5 [-], while the opposite side is submitted to 
interior loading conditions 6'int = 24 [°C] and ipi^t = 0.8 [-]. The solution of 
the time-dependent problem in Eq. ( 129|) also involves a discretisation of the 
time domain T into T = 151 time steps and hence the PCE-based surrogate 
model consist of x T = 12, 080 PCEs for the temperature, and the same 
for the moisture. 

In order to describe the accuracy of such a surrogate model, let us define 
the MC estimate of the error expectation e{u) as a relative difference between 
two response fields tt"^ and over the discretised spatial and time domain 

as 

V„=i t=i / 

The quality of a PC-based surrogate model depends on the number M 
of eigenmodes involved in KLE describing the fields of material properties 



13 



as well as on the degree of polynomials P used in the expansion Eq. 
Figure [2] shows the error estimate e{u) computed for different numbers of 
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Figure 2: Errors in approximation of the temperature (a) and the moisture (b) field 
induced by PCE and KLE as a function of number of eigenmodes. 

eigenmodes M and for the polynomial order P = 2 and P = 3. Here, the 
response fields it" are computed by the FEM based on one realization of the 
KLE of the parameter fields (further shortly called FE simulations) and the 
response fields are obtained by evaluation of the constructed PCE in the 
same sample point. In order to distinguish the portion of error induced by 
the KL approximation of the parameter fields, the estimate e{u) is computed 
once for the FE simulations using all M = 120 (dashed lines), and once for 
the FE simulation using the same number of eigenmodes as in the constructed 
PCE (solid lines). In other words, the solid lines represent the error induced 
by PC approximation and the difference between the solid and corresponding 
dashed line quantifies the error induced by the KL approximation of the 
parameter fields. 

Figure [3] represents the same errors e{u) as Fig. |21 but this time with 
respect to the computational effort needed for the computation of PC coeffi- 
cients. Regarding the obtained results, we focus our following computations 
on the KL approximation of the material parameters including M = 7 eigen- 
modes and a PCE of order P = 2 providing, at reasonable time, sufficiently 
good approximation of the model response, namely of the temperature field 
where the errors are more significant. 



^We assume the full PC expansion, where number of polynomials R is fully determined 
by the degree of polynomials P and number of eigenmodes M according to the well-known 
relation R = (M + P)!/(M!F!). 
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Figure 3: Errors in approximation of the temperature (a) and the moisture (b) field 
induced by PCE and KLE as a function of computational time needed for a PCE con- 
struction. 



For a more detailed presentation of the PCE accuracy, Fig. H] compares 
the model response in one node of FE mesh (the node No. 1 at Fig. [5]) at the 
time t = 400 [h] obtained by the FE simulation and by the PCE as a function 
of the first stochastic variable ^i. 




(a) (b) 

Figure 4: Detailed comparison of the temperature (a) and moisture (b) with their PC 
approximation as functions of the first stochastic variable . 



6. Bayesian updating procedure 

In the Bayesian approach to parameter identification, we assume three 
sources of information and uncertainties which should be taken into account. 
The first one is the prior knowledge about the model/material parameters 
qi{uj) defining the prior density functions. In our particular case, we know 
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that all the identified parameters are positive-definite and the log-normal 
random fields with the statistical moments given in Tab. [T] are suitable for 
a description of the prior information. In fact, they are maximum entropy 
distributions for this case. We describe the material parameters using the 
KLE which is fully defined by a finite set of standard Gaussian variables 
^ = [^i(w), . . . ,C,m{^)] with the probability density function (pdf) p^{$) and 
thus, the updating procedure can be performed in terms of ^ turning them 
into non-Gaussian variables. 

Other source of information comes from measurements, which are vio- 
lated by uncertain experimental errors t{u}). Last uncertainty Cj arises from 
imperfection of the numerical model, when for example the description of 
a real system does not include all important phenomena. However, it is a 
common situation that the imperfection of the system description cannot be 
distinguished from measurement error e and the modelling uncertainties Co 
can be hidden inside the measuring error e{u). Then we can define the pdf 
Pz{z) for noisy measurements z{u}). 

Bayesian update is based on the idea of Bayes' rule defined for prob- 
abilities. Definition of Bayes' rule for continuous distribution is, however, 
more problematic and hence jsi Chapter 1.5] derived the posterior state of 
information 7r(^, z) as a conjunction of all information at hand 

7r{tz) = f,p^{^)p,iz)piz\^), (31) 

where k is a normalization constant. 

The posterior state of information defined in the space of model param- 
eters $, is given by the marginal pdf 

7r^(^) = En z)] = Kp^U) [ Piz\^)p.iz) ndcu) = np^i^) L(^), (32) 

where J? is a set of random elementary events u and measured data z enters 
through the likelihood function L{^), which gives a measure of how good 
a numerical model is in explaining the data z. 

The most general way of extracting the information from the posterior 
density vr^(^) is based on sampling procedure governed by MCMC method. 
For more details about this approach to Bayesian updating of uncertainty 



in description of couple heat and moisture transport we refer to [11|. In 
this paper, we focus on the comparison of the posterior information obtained 
from the sampling procedure using directly the computationally exhaustive 
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numerical model f ll6p on one hand and using the PC approximation of the 
model ( 12 5 p on the other hand. 



7. Numerical results for the Bayesian update 



Due to the lack of experimental data, we prepared a virtual experiment 
using a FE simulation based on parameter fields obtained by the KLE with 
7 eigenmodes so as to avoid the error induced by KLE, which is mainly the 
subject of the work presented in ll|. A related set of random variables ^ 



is drawn randomly from the prior distribution and stored for a purpose of 
latter comparison with the prior and the posterior state of knowledge. The 
resulting temperature and moisture fields considered as a so-called "true 
state" or simply the "truth" are shown in Fig. \5\ According to the 
values of temperature and moisture are measured in 14 points (see Figs. 
(a) and (c)), and at three distinct times (see Figs. |5](b) and (d)). Hence, the 
observations z consist of 84 values. To keep the presentation of the different 
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Figure 5: Virtual observations: (a) and (c) spatial arrangement of probes; (b) and (d) 
temporal organization of measurements 



numerical aspects of the presented methods clear and transparent, we focus 
here on a quite common and simple case, where modelling-uncertainties are 
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neglected and measurement errors are assumed to be Gaussian. Then the 
hkehhood function takes the form 

m = '^exp (-i (Y(^) - zfC-l - z)^ , (33) 

where Y{^) is an observation operator mapping the model response u given 
parameters ^ and loading / to observed quantities z. Cobs is a covariance 
matrix representing the uncertainty in experimental error, which is obtained 
by perturbing the virtual observations by Gaussian noise with standard de- 
viation for temperature ag = 0.2 [°C] and for moisture = 0.02 [— ] so 
as to get 100 virtual as an input for the covariance matrix evaluation. In 
order to be able to compare the posterior state with the true state also in 
terms of model parameters ^, we assume an artificial situation where the 
observed quantities z correspond exactly to the true state of temperature 
and moisture. 

The Bayesian update was performed using Metropolis-Hasting algorithm 
and 100, 000 samples were generated in order to sample the posterior density 
f PT]) over the variables ^ = (^i ■ ■ ■ ^m=7)- The truth state, prior and posterior 
pdfs obtained by the FE simulations and using the PCE are plotted in Fig. |6l 
One can see that the error induced by PC surrogate of model response are 



fPosterior - FE(M = 7) 




(a) (b) 

Figure 6: Comparison of pdfs (a) for the separate variable and (b) for the pairs of 
variables. 

negligible in terms of the resulting posterior densities. Figure [6] also demon- 
strates the fact that the variables ^ being a priori standard Gaussian should 
not be a posteriori Gaussian. 

During the sampling procedure, we stored also the corresponding values 
of parameter fields and response fields in order to obtain their posterior state 
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of information. As a result, Fig. [7] shows the comparison of the truth, and 
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Figure 7: Comparison of pdfs for material properties in FE node 7: (a) the thermal 
conductivity of dry material Aq and (b) water vapour diffusion resistance factor fj,. 

prior and posterior pdfs for two material parameters Aq and /i in the top-right 
corner FE element, and similarly. Fig. [8] presents pdfs for the temperature 
and moisture in FE node 7 at 400[h] (i.e. at the 151-th time step). 
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Figure 8: comparison of pdfs (a) for the temperature and (b) for the moisture in FE node 
7 at 151-th time step. 

We should note that the similarity of the prior and the posterior pdfs 
for moisture in Fig. |8] is probably caused by the very slight influence of the 
studied material parameters to the moisture value or more precisely, the prior 
standard deviations were very small. 

Beside the comparison of the PCE accuracy, we also compared the time 
necessary to generate the samples. In case of PCE, the total time also in- 
cludes the time of PC coefficients computation. Particular comparison of 
computational time needed by FE simulations and by PCE evaluations is 
demonstrated in Fig. |9l 
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Figure 9: Comparison of time necessary for evaluation of samples. 
8. Conclusions 

The presented paper presents an efficient approacli to propagation and up- 
dating of uncertainties in description of coupled heat and moisture transport 
in heterogeneous material. In particular, we employed the Kiinzel's model, 
which is sufficiently robust to describe real-world materials, but which is 
also highly nonlinear, time-dependent and is defined by 8 material parame- 
ters difficult to be estimated from measurements. The updating procedure 
starts with the prior information about the parameters' properties such as 
positive-definitness and second order statistics. Heterogeneity of the material 
under the study is taken into account by describing the material properties 
by random fields, which are for a simplicity considered as fully correlated. 
Then, the corresponding correlation lengths are assume to be known as an- 
other a priori information. In order to limit the number of random variables 
necessary to describe the material, the random fields are approximated by 
Karhunen-Loeve expansion and hence, all the remaining uncertainties are 
described by a set of standard Gaussian variables whose number is given by 
the number of eigenmodes involved in KLE. 

These uncertainties are then propagated through the numerical model 
so as to provide a probabilistic characterization of the model response, here 
the moisture and temperature fields. Simultaneously, the other information 
including uncertainties coming from the experimental measurements is used 
to update the prior uncertainties in the model parameters. In order to imitate 
the experimental measurements, a virtual experiment is prepared together 
with the relating uncertainties given by a covariance matrix. The Markov 
Chain Monte Carlo method is then employed so as to sample the posterior 
state of information. 
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The primary objective of the presented paper is to accelerate the sam- 
phng procedure. To this goal, a polynomial chaos-based approximation of the 
model response is constructed in order to replace computationally expensive 
FE simulations by fast evaluations of the PCE during the sampling. In par- 
ticular, the PC coefficients are obtained by an intrusive stochastic Galerkin 
method. It is shown that the resulting approximations exhibit high accuracy 
and the related posterior probability density functions are sufficiently precise 
as well. Finally, the comparison of the computational effort confirmed the 
large savings in case of PC evaluations. 

While the acceleration obtained by the presented procedure is significant, 
it can be still unfeasible for very large problems. Our future work will be 
focused on the elimination of the MCMC sampling procedure itself by the 
update directly in terms of parameters of probability density functions as 
proposed in [36| . 
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